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Abstract 

The Gaudin integral equation for the ground state of a one-dimensional ^-function attractive 
spin- 1/2 fermions is solved in the form of power series. The first few terms of the asymptotic 
expansions for both strong and weak coupling cases are calculated analytically. The physical 
quantities such as the ground state energy are expressed in terms of a single dimensionless 
parameter 7 = c/D, where c is the coupling constant and D is the number density. The re- 
sults agree with those obtained from the perturbation calculations, which include the one in 
the classical electrostatics originally by Kirchhoff. In the strong coupling limit, the connection 
to the solutions of the Licb-Linigcr integral equation for the ground state of a one- dimensional 
(5-function repulsive bosc gas is shown explicitly. 

KEYWORDS: (5-function gas, one-dimensional integrable systems, BCS-BEC crossover, 
Lieb-Liniger integral equation, Gaudin integral equation, Tonks-Girardeau gas, hard-core bose gas 

1 Introduction 

We discuss a one-dimensional system of quantum N particles with ^-function interaction in a 
periodic box of length L. The Hamiltonian of the system is 

N . 2 a2 

* = -E2^+«E^-*i). (i-D 

3=1 3 

where m is the mass and c is the coupling constant. Hereafter, we choose a unit, h = 1 and 
2m = 1 to simplify the expressions. 

The system I|1.1J) has been studied extensively by use of the Bet he ansatz method. 1-9 For 
N repulsive (i.e. c > 0) bosons, the ground state is described by the integral equation for the 
distribution function p{k) of the quasi-momenta, 1 

2tt tt J_ k c 2 + (k -q) 2 ' y ' 

where A' is a cut-off momentum, which is fixed by the normalization condition, 

p(k)dk = D. (1.3) 

K 



*iida@monet.phys.s.u-tokyo.ac.jp 
^wadati@monet.phys.s.u-tokyo. ac.jp 



1 



Here, D = N/L is the number density of the particles. The interval [— K, K] is filled with 
quasi- momenta k. The ground state energy is given by 



E B 



L 



k z p{k)dk. (1.4) 

K 



We call ((1.2(1 the Lieb-Liniger integral equation. 

In 1967, Yang 4 proposed a method, which we now call the generalized Bethe ansatz, or 
"nested" Bethe ansatz, to study the system ()l.l|l with arbitrary internal degree of freedom. For 
N attractive (i.e. c < 0) spin-1/2 fermions, the ground state with total spin is described by 3 

n J-q c +{q- q ) 

Here, f(q) is the distribution of spins with spin- wave rapidity g, and a cut-off rapidity Q is fixed 
by the normalization condition for f(q): 



f{q)dq = D, (1.6) 



while the ground state energy is given by 



f f- C *-)f{q)dq. (1.7) 



L 7T 7_ Q V 4 
For our purpose, it is convenient to introduce the effective ground state energy, 

^f^^ + ^D= 1 -[ Q qlf(q) d q. (1.8) 

L L 4 7T J_ Q 

We call (ll.5fl the Gaudin integral equation. 

By use of a method proposed in the previous paper 10 to solve 1(1.2(1 . we explicitly calculate 
the first few terms of the asymptotic expansions of /(g) for both strong and weak coupling cases. 
It is essential to use the single dimcnsionlcss coupling constant 7, defined by 

c 2mc 

1= d = Wd- (L9) 

Combining the results with those in the previous papers, 10, 11 we discuss in particular four cases, 
7 — > —00, 7 — > — 0, 7 — ► +0, and 7 — ► +00. We also study the physical quantities, such as the 
ground state energy and the chemical potential, as functions of 7. 

This paper is organized as follows. In § we solve (jl.5(l for large and small | — ^ | . The 
similarities and the connections of the obtained expressions in the case of 7 — > —00 for p. 5(1 and 
7 — * +00 for H.2fl are investigated in § O The recent discussion on "BCS-BEC crossover" 12-14 
is made clear for these integrable models in one-dimension. For the small I7I case, it is known 
to be difficult to solve 1(1.2(1 and 1(1.5(1 . Analogies with the perturbation calculation in the 
classical electrostatics originally by Kirchhoff are discussed in § 0] The last section is devoted 
to the conclusions. To concentrate on physical discussions, most of the calculations are given in 
Appendices 



2 Solutions of the Gaudin Integral Equation 
2.1 Formulation 

Following a method presented in the previous papers, 10 ' 11 we solve the Gaudin integral equa- 
tion 1(1.5(1 . In terms of new variables, 

q = Qx, c=-QX, q' = Qy, (2.1) 

f(Qx)=F(x), (2.2) 
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Ij 1.5(1 becomes 



a r 1 
F(y) + - \ 
n J-i 



F(x) dx 



= 2. 



(2.3) 



A 2 + (x - y) 2 

We also refer to ()2.3(l as the Gaudin integral equation. As mentioned, we use the number density 
D and the dimensionless coupling constant 7, 



D = 



N 
c 

1= D- 



(2.4) 
(2.5) 



The large (small) A corresponds to the large (small) I7I 
By expanding 



F(x) = J2 



2n 
(l n X , 



(2.6) 



we have 



F(x) 



1 dX X 2 + {x _ y) 2 
00 00 P n n2m+l 



EEE 



n— p— m— 



(2m)l (p — m) 



Qm+p 
Q s m+p 



m+n ( ._ 1 y+i A 2(;-l) 

^ 2(m + n - /) + 1 



(_l)m+»A 2 ( m +")- 1 arctan- 

A 



, 2 '' 



(2.7) 



/here s = A 2 . Substituting l|2.6|) and 1(2. 7jl into l|2.3|) . we get 



a vV + - 2^ 2^ 2^ ( 2m )!(p-m)! 



p=0 



Qm+p 

g s m+ P 



E 



(-ly+iA 2 ^- 1 ) 



^— ' 2(to + n — 
. 1=1 y 



+ (_!)«+" A 2(w+n)-l arctan 

O + i 



, 2 P 



(2.8) 



In order that (|2.8|) holds for any y, — 1 < y < 1, we equate the terms of the same power of y 2p . 
At y° and y 2p (p > 1), we respectively obtain 



2 1 
do H — arctan — • ao 



7T 

2A 

7T 



E M E ^ n 1 ; + (-i)'^ 2 "- 1 



n=l 



a - 2(n-Z) + l 

a„ • 2 2m d m+p 



arctan ■ 



p 7T ^ ^ (2m)! (p - m)! ds m +P 

n=p+l m=0 v y v ^ ' 



m+n 

E 

. i=i 



(-l) 1+1 s 



2(m + n - /) + 1 



00 p 



+ vEE 



a„ • 2 



2i» 



Qm+p 



it ^ (2m)\ (p - m)\ ds m +P 

n=0m=0 v y ^ ; 



(-l) m + n A 2(m+n)_1 arctan 



A 



0. 



(2.9a) 



(2.9b) 



In this way, a solution of ((2.9() gives a solution -F(x) of (|2.3(1 . 

Until now, all the expressions are exact and no approximation is employed. We summarize 
here the method of solution: 

1. Solve (|2.9() to determine {a n } of (|2.6|l in the case of large A and small A. 
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2. Use 1(1.61) with 1)2. which can be rewritten as 

ttA 

7 = 



(2.10) 



to obtain A as a function of 7. 

3. Calculate other physical quantities from the expressions of (|l.(j|) . 12.1(1 . (|1.8|) and thermo- 
dynamic identities 2 in the following form; 



(2.11) 
(2.12) 
(2.13) 
(2.14) 



D 


=17 






Ef 




L 




/' 


= D 2 ( 



■ 2 F(x) dx, 

7 d 7 /' 

de 1 9 d 2 e x 
" 2 ^ + 2^V 



where e = E-p/(ND 2 ), /i and 1; are the scaled energy of the ground state, the chemical 
potential and the sound velocity, respectively. 



2.2 Strong coupling case 

We consider the large A case. In this case, only a few a„'s are enough to have good approxima- 
tions. Indeed, it can be shown that a n = 0(1/A 2n+1 ) for n > 1. Therefore, we only keep ao, a±, 
and set a n = (n > 2), 

F(x) = ao + aix 2 . 



Equations 1)2.9(1 give 



2 1 2/11 

a + - arctan - ■ a + — - - - — 
7r A 7rA V 3 5A^ 



ai = 2, 



2 2 



and consequently, 



a = 2 1 - — 



2tt 2 - 24 



ttA tt 2 A 2 3tt 3 A 3 



ai 



ttA 3 ' 



Using 1)2.15(1 with 1(2.17(1 in 1(2.10(1 . we obtain A as a function of 7. The result is 



A 



2(27+1) 



The expressions of ao and a\ in terms of 7 are 

2 7 +l 



a 



a 1 



7 



1 - 



6(2 7 + l) 3 



2 

7T Z 7 

6(27 + iy 



'2(2 7 + 1)3' 



(2.15) 

(2.16a) 
(2.16b) 

(2.17a) 
(2.17b) 

(2.18) 

(2.19a) 
(2.19b) 
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The ground state energy is calculated from (|2.12|) . 



l -^{WTi) [ 1+ wTTW*)- (2 - 20) 

The cut-off rapidity Q, the chemical potential \i and the sound velocity v are given by 
ttD 7 



2 2 7 + l V 1 + 6(27+l)V' (2 ' 21) 
S-**P^(^)'( 1+ »*J\, (2.22) 



12 27 + fV 2 7 + 1 / V 5(2 7 + f) 2 (67 + f) 
,27TT) ( 1+ 3(27+1)0' 



^Vfl + ^T^V (2-23) 



where we have introduced the effective chemical potential \x , = /i + (1/4)D 2 7 2 . 
2.3 Weak coupling case 

This is the small A case. The lowest order approximation of 12.9fl is 

2a = 2, (2.24a) 
2a p = 0, p > 1 . (2.24b) 

A derivation of l|2.24(l from ({2.9(1 is the same as shown before. 10,11 Another derivation is given 
in Appendix where we treat the integral in a more direct manner in the small A case. 
Equations l|2.24|l simply gives 

F(x) = 1. (2.25) 

We list here the analysis based on 1)2.25(1 : 

A = -- 7 , Q=^, (2-26) 

poff _2 _2 

f^=D 3 — , // eff = D 2 — , _ j^jj (2.27) 

L 12 M 4 ' v y 

It should be noted that, in contrast to the solution of the Lieb-Liniger integral equation 11.2fl . 
all the coefficients a n are finite at A — > 0. 10 This nice analytic property can be seen from 1(2.3(1 . 
since the kernel A/{7r[A 2 + (x — y) 2 )]} becomes an expression of 5(x — y) in the limit A — *• 0. 
The second order approximation of ((2.9(1 is 

2a p + -A V ^j" 1 a n = 25 p0 , P > 0. (2.28) 
7r ' 2n — 2n — 1 



n=0 



Here, 8 p o is the Kronecker's delta. A derivation of l|2.28() is also given in the previous papers 10 ' 11 
and in Appendix lAl We set 

a n = S n0 + a^> + o(A), n > 0, (2.29) 

where 8„o stands for the solution of (|2.24|) and denote the corrections of the A 1 order. 
Substitution of (|2~2^|) into (|2~2"5|) yields 

2aW - -A = 0, p > 0, (2.30) 

7T 



which leads to 

F(x) = 1 + - Y x 2n + o(A) 



A °° 



TT 

n=l 



^ + (A). (2.31) 
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This agrees with the approximation shown in the Gaudin's book. 7 From (j2.10|) . F(x) is inte- 
grated to give 



2b,^ 1 



7 = -7rA^2 + — A J , (2.32) 
where the divergent sum b is defined by 

oo 

b = Y- . (2.33) 

^ 2n + 1 ( ' 

n=0 

Equation (|2.32|) is solved to obtain 7 as a function of A: 

A = --7(1-^7). (2-34) 
We proceed to calculate the physical quantities. From 112.11(1 and (j2.34(l . we have 



ttD ( 2b 

~2~ { ~ 7T 



1 - —7 ) . (2.35) 



The ground state energy is obtained from l|2.12|l , (|2.31|) and (|2.35|l : 



Ef _ Q 3 (l 2A^, 1 

\ n=0 



,2 



12 



12 ^ + ^ 7 £j(S7Ti)(2„ + 3)y 

+ ;«,). (2.36) 

This recovers the previous result 16 where the integrals are directly treated without obtaining 
the explicit form of F(x). It is remarkable that we observe the cancellation of the divergent 
sums as before. 10,11 That is, the physical quantities are expressed in power series of 7 without 
any divergent sum. The chemical potential and the sound velocity can be calculated from the 
formulae ({2~T5j) and (|2~T4"|) . 

^ = D^(l + -i 7 ), (2.37) 



v = irD ( 1 + ^ ) . (2.38) 



7T- 



The limit 7 — * —0 gives the well-known results for free fermions. 

In fig. E w e plot the solutions Ij2.15jl with Ij2.19jl and (|2.31|) with (|2.34|) , for several values 
of 7. For comparison, we also plot the numerical solutions based on the classical quadrature 
method. 1,17 In this figure, we choose a unit, D = 1 so that Q/D takes a value between ir/4 
and 7r/2 (see (|2.21|l and H2.35|l '). For a fixed Q and — 1 < x < 1, q and /(g) are determined 
by (|2.1(l and (|2.2(l . respectively. Numerical values of 7 are obtained through the relation 7 = 
—\Q/D. In the weak coupling case, it is fair to say that the error is large near the end-points, 
\q\/D ;$ Q/D. This can be observed analytically when the solution H2.31(l is substituted into 
the integral equation (|2.3|) . 



3 The BCS-BEC Crossover 

Recently, there is a renewed interest in the crossover from a Bardeen-Cooper-Schrieffcr (BCS) 
superfluid with Cooper pairs to a Bose-Einstein Condensate (BEC) of molecules composed of 
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Figure 1: The normalized distribution function of spins f(q) in the ground state for various values of 
7. Numerical results are indicated by points (+, x etc.), while the analytical ones, ()2.15[) and l|2.31|) . 

are by lines ( , etc). Numerical solutions are obtained by a classical method. 1 ' 17 For small 

7, analytical solutions are shown for the intervals [-Q/D + |7|/4, Q/D — |7|/4] of q/D, which 
correspond to [— 1 + A/4, 1 — A/4] of x. For intermediate 7, 7 = —1.136 and 7 = —2.239, only the 
numerical results are shown. For large 7, analytical and numerical results are indistinguishable. 
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tightly bound fcrmion pairs. 13, 18 The experimental realizations of low-dimensional trapped bosc 
gases at ultracold temperatures also inspired the controversy even in one-dimension. 12 ' 14 The 
analytic solutions in § [21 of the Gaudin integral equation and in the previous paper 10 of the 
Lieb-Liniger (LL) integral equation enable us to make the quantitative discussion. To avoid 
confusion, physical quantities, such as the mass m, the number of particles N, the number 
density -D, the coupling constants A, 7, and the ground state energy E are denoted with the 
subscripts F and B for a fermi gas with spin-1/2 and a bose gas, respectively. 
The dimcnsionlcss form of the LL integral equation i|1.2|l reads 1,10 



9(y) - ^ f 
k J-i 



g{x)dx 1 



\l + {x-yf 2tt' 
where 

g(x) = p(Kx), k = Kx, 2m B c = K\ B . (3.2) 



The relation between Ab and 7b is 



2to b c A b , . 

7B = ~^b" = T 1 ' ( } 

g(x) dx 



where D B = N-q/L. The ground state energy is given by 

5= <"> 

As already mentioned in the Gaudin's paper, 3 strong coupling fcrmions described by the Gaudin 
integral equation 12.3fl may be regarded as the repulsive bosons of dimers. To examine whether 
it is indeed true, we establish a relation between the exact solutions of the Gaudin integral 
equation (|2.3() and the LL integral equation 1|3.1|) . In the following, we concentrate ourselves 
on the strong coupling case, i.e. the large Af case for fermions and large Ab case for bosons 
are considered. For explicit comparison, we summarize here the results shown in § 12.21 and the 
previous papers. 10,11 

For — 7f 1 (fcrmions with strong attraction), 



2 2 7F + 1 V 6(2 7F + l) 3 
Ef _ Dl n 2 ( 7F \ 2 ( 4tt 2 



L 2m F 12\2j F + lJ V 1 + !5(27f + I) 3 ) ' ^' 5 ^ 

ff Dj , 2 6 7F + 1 / 7F V ( Sn 2 \ 

fiF 2m F 12 2 7F + 1 ^2 7F + iy ^^ + 1)^ + 1)/' [ ' 

2to f V 2 7f + 1/ V 3(2 7F + l) 3 y 1 j 
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For 7b 3> 1 (bosons with strong repulsion). 



7B + 2 / 2if_ 

27T 7 B I, 3 7B (7B + 2)V ^ B D i 



= J K=r~ ( 1 - 57 ~. , ^ ) - ZT^fc 2 - 0.6a) 



AB = ^fl-^^), 0.6b) 



7T 



3( 7 B + 2) 3 



7 B + 2 V 3( 7B + 2) 
E B _ D% t: 2 f 7 B V A , 32tt 2 



L 2m B 3 V7B + 2; V 15( 7B + 2) 3 

D B 3 7B + 2 7T 2 ^ 7 B \ 2 f, , 64tt 2 



1 + 777 ' ( 3 - 6d ) 



MB 2m B 7b + 2 3 V7B + 2J ( v 1+ 5(7b + 2) 2 (3 7 b + 2) / ) ' 
27rD B / 7b Y fi 16tt 2 



Vb= *^{^T2) { 1 + 3(^TW)- (3 ' 6f) 

Identifying Afp fermions with mass m-p with N B bosons of dimers with mass m B , we have 

N F D F 
m B = 2m F , AT B = — , D B = — . (3.7a) 



Hence, from ljl.9|l we have 



Substitution of (|3.7|) into (|3.5dJ) gives 



2 ' 2 



7b=47f. (3.7b) 



£f ff (2£ b ) 3 tt 2 / 7b/4 V /, , 4 ^ 2 



L 2m B /2 12 Vtb/2 + 1/ V 15(tb/2 + l) 3 
Dl 7T 2 f 7B V ^ , 32tt 2 \ E B 



2m B 3 V7B + 2y V 15( 7 b + 2) 3 ; L 



(3.8) 



The effective ground state energies are exactly the same through the relation l|3.7|) in the strong 
coupling case. We can easily verify that other physical quantities are also smoothly connected 
through the relation H3.7J) at the point 1/tf = 1/7b = 0- In fact, in the case being considered, 
equations Q2.3|l and can be treated in a unified way, if we regard A as a parameter and 

D as a normalization condition for these integral equations. We perform such calculations in 
Appendix IB1 

We note that the quasi-momcnta for bosons and fermions are related as 

K = 2Q, k = 2q. (3.9) 

These relations are consistent with the assumptions that lead us to the Gaudin integral equation. 
For negative c, quasi-momentum {k} of the fermions (not to be confused with k for the repulsive 
bosons in l|3.9|) ) become imaginary. In the thermodynamic limit, two k s with a real part g, 
which arc complex conjugates of each other, form a bound state of two fermions with the 
binding energy c 2 /2 (see (|1 . 7|) 'l . 9 This kind of argument usually starts with the c = —0 limit. 
Remarkably, it holds also in the strong coupling case. 



Equations (|3.5d|) and l)3.6d(l need some explanation. Since we are seeking for the solution 
of 0(1/A 3 ), the expansion must be in powers of either 1/A ~ 1/(7b + 2) ~ 1/(tf + 1/2) or in 
I/7B — 1/7f- However, it occurs that the solution of the first term should be (tb/(tb + 2)) 2 — 
( I Tf |/(| Tf I ^ 1/2)) 2 , not be simply 1. This factor can be counted as the reminiscent of the 
effective "diameter" a of each particle, in the limit of the hard-core bose gas, or, in other words, 
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the quantized version of the classical Tonks gas. 15 In fact, the energy at zero temperature for 
the hard-core bose gas is calculated as 19 



E hc 7T 2 / 1 



L 2rriB 3 \ 1 — oDb 
Besides, equating the scattering phase shifts in the low energy limit, we find 2 

12 11 



2ms c 2mp c 



(3.10) 



(3.11) 



With the relation (|3.11(l . we see that Q3.1()|l and the first terms of l|3.6d|) and l|3.5d|) become 
exactly the same. The positiveness (negativeness) of the diameter a for the attractive fermions 
(repulsive bosons) is due to its statistical property. The negativeness for the repulsive bosons 
can also be interpreted as the transparency effect for a finite c > 0. It is easy to verify that the 
chemical potential /ihc and the sound velocity Uhc of the hard-core bose gas 19 are obtained by 
the substitution of (|3.11() with i|3.7|l into cither i|3.5c|l and (I3.5f|) or (|3.6cfl and i|3.6i|) . Discarding 
the second term in the last brackets, we obtain 

P 2m B 3 1 - aD B \ 1 - aD B J 
2tt£> b ( 1 V 

^ = ^\Y^D B ) ' (3 ' 13) 

Note that /ihc corresponds to 2/zf?. 

Now we are ready to describe the BCS-BEC crossover in one-dimension. The discussion 
above for smooth connections in the strong coupling cases justifies a link between the two 
systems. In view of the physical quantities, we can map two fermions of opposite spins forming 
a bound state whose quasi-momenta are g± ic/2 onto one boson with quasi-momcntum k = 2q. 
In such a way, the effective ground state energy per volume, the chemical potential and the 
velocity of sound can be identified. If we take 7f as the coupling constant common to the both 
systems and vary 1/7f from — 00 to +00 through 0, we find the following asymptotic expressions 
of the effective energy: 



(3.14) 




where cff is the Fermi energy of the free fermions. 

1 / irDp \ 



6ff 



2m F V 2 



I . (3.15) 



The meanings of the above three expressions are clear. They are the energy of the Hartree- 
Fock approximation for weakly interacting fermions, 16 the energy of either the hard-core bose 
gas or the spinless non-interacting fermions, 26 and the energy due to the Bogoliubov theory 
for weakly interacting bosons, 1 for — I/7F 3> 1, I/7F ~ 0, and I/tf 3> 1, respectively. The 
asymptotic solutions of the LL and the Gaudin integral equation reproduce the results of these 
three "phases" rigorously. Two crossovers for the intermediate value of 7f are by no means a 
phase transition. The solutions are continuous and analytic functions of 1/tf except the two 
points —00 and +00. Those can be proved mathematically from the properties of the integral 
equations. 
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Figure 2: Asymptotics of the ground state energy in units of the Fermi energy of the free 
fermions as a functions of the inverse of the coupling constant, 1/jf- 
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Figure 3: Three regimes of one-dimensional quantum gas: for 7f < 0, spin-1/2 fermi gas and for 
7f = 7b/4 > 0, 5- function bose gas. 
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Let us explain these three regimes in detail. When — 1/tf 3> 1, the ground state is a BCS- 
like state with Cooper pairs. 16 The effective size of this bound state is much larger than the 
average interparticle distance. The existence of the thermodynamic limit is guaranteed by the 
Pauli exclusion principle; i.e. there is no bound state of more than two fermions. 

The intermediate regime of the hard-core bose gas around 1/tf = should be distinguished 
from the Girardeau gas, which corresponds only to the limit I/7F = 0. Indeed, when we 
compare our result with the numerical one, as in fig. it traces well the sudden change around 
1/Tf = 0. In this region of the bosonic side, the strong mutual repulsion compels the bosonic 
particles to exhibit the fermionic properties due to one-dimensional kinematics. On the other 
hand, even the infinitely strong attraction cannot causes the fermions to condense. Namely, the 
distribution and the energy spectrum of the tightly bound fermion pairs can be identified with 
those of spinlcss hard-core bosons with a positive "diameter" . 

In the regime of weakly interacting bose gas, the lowest order solution of the distribution 
p(k) is 1 ' 10 



P (k) 



1 



(DhF 



.1/2 



(3.16) 



27tDf7f 

This indicates that the weaker the interaction, the narrower the quasi-momentum distribution 
region and the higher the peak of the distribution. In the recent experiment this behavior was 
observed for 73 = 0.5. 21 In the limit I/7F — ► +00, all the particles condense into the zero 
quasi-momentum state. 

In the same manner as we have done in the effective ground state energy, we list here the 
asymptotic expressions of the velocity of sound and the effective chemical potential: 



1 + 2* 



VFF 



7f 



2 7F 
1 1/2 



1 



2tt 2 



1, 

1 1 1/2 



3(2 7F + If 



-1/7F > 1, 
1/7F ~ 0, 
1/Tf > 1. 



(3.17) 



where vff is the Fermi velocity of the free fermions, 

7tD f 
2mp 



VFF 



(3.18) 



and 



£ff 



4 

1 + 

1 67F + 1 



7f 



3 2 7F + 1 \27F + 1 
2 / 2 
7f 1 



8tt 2 



5(2 7f + 1) 2 (6 7 f + 1) 
These expressions are plotted in fig. 01 and fig. as functions of I/7F. 



1/2 
7f 

k 7T V 7T 



-1/Tf > 1, 
1/7F ~ 0, 
1/7F > 1. 



(3.19) 



4 Analogy with the Circular Plate Condenser 

An intriguing fact is that the Gaudin integral equation and the LL integral equation also appear 
in a completely different subject, the circular plate condenser. 5, 7,22 We consider the condenser 
of the radius 1 whose two plates are separated by a distance A. Cylindrical coordinates are 
taken as follows: the radial distance from the common axis of the disks is taken to be p, the 
distances from the disks to be £, respectively. E. R. Love 22 (1949) proved that the potential 
4>(p, C, C) due to the disks is expressed as 

<KM,C') = — / I m ± " un\h(t)dt, (4.1) 

" J-i\[p 2 + (( + it) 2 } 1/2 [p 2 + (C - it) 2 ] 11 2 J 
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Figure 4: Asymptotics of the sound velocity v-p in units of the Fermi velocity of the free fermions 
uff as a function of the inverse of the coupling constant, 1/tf- 
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Figure 5: Asymptotics of the chemical potential /zp ff in units of the Fermi energy of the free fermions 
£ff as a function of the inverse of the coupling constant, 1/tf- 
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where h(t) is determined by 



A r 1 
h(t)±- \ 
k J-i 



h(s) ds 



(t - s) 2 + A 2 
The capacitance C of the condenser is given by 

1 



= 1. 



C = 



2tt 



h(t) dt. 



(4.2) 



(4.3) 



In the above expressions, the upper signs refer to the case of equally charged disks at potential 
Vq, and the lower signs to that of oppositely charged disks, namely, disks at potential Vq and 
— Vq. The Love equation (|4.2|l is exactly in the same form as the Gaudin and the LL integral 
equation: to be precise, h(x) in (|4.2(l corresponds to (l/2)F(x) and 2irg(x) in (|2.3|) and JsQJj 
respectively. 

In the case of oppositely charged disks (lower signs), the solution of l|4.2[l is calculated for 
small A (small separation). By the heuristic argument, 23,24 away from the edges hit) to the 
o(l) order is written as 



^) = ^(l-* 2 ) 1/a + (1 "^" 1/2 



, 167re , 1 + t 
log— -Hog — 



and with (|4.3|l . the capacitance is 



1 1 16tt 



This result is originally by G. Kirchhoff (1877). 25 
On the other hand, we have 10,20 



and 



C = 



2ttA 



1 oo oo 



(2m + 1)!! (2m - 1)!! (2n-l)!! 
[(2m)!!] 2 (2n)!! (2n-2m-l) ' 



/ q(x) dx = — — + -b, where b = 

7 i 4A 27r 4 ^ 



(2m- 1)!! 



(2m) 



We regularize the divergent sum b as follows. We set 



m = E 



(2m - 1)!! 
(2m)!! 



,2m 



-K{k). 
7T 



(4.4) 



(4.5) 



(4.6) 



(4.7) 



(4.8) 



The function A"(fc) is the complete elliptic integral of the first kind with the modulus k. Asymp- 
totic behavior of K(k) is 

as k 2 / 1. 



If we choose 1 — fc 2 = A/(7re) in 14.9|l with (|4.8|) . we get 

b = lim b(k) = — log — . 

fc2/l 7T A 

Using the expression (|4.1()() . we may rewrite 1)4. 6JI as 



1 



1 OO OO 

2=0 m=0 



(2m- 1)!! 



(2m)!! 



OO OO 



(2/ - 1)!! 2Z 



4tt 



i=0 m=0 



(20! 



(2m - 1) 



(2m) 



(20!! 

1 ^ 



2m- (2Z- 1) 



— (1- x 2 ) 1/2 + — 
2nX [ ' + 4^ 2 



, 167re , 1 + x 
log — x log 



1 



l-a?y (l-a; 2 ) 1 / 2 



(4.9) 



(4.10) 



(4.11) 
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Here, we have used the power series formulae, 

E^^*^ 1 -* 2 ) -1 ' 9 ' N<1 ' (4 - 12) 

p=o \ P>-- 

and the following summation formula, 

1 2 (2QH y> 1 (2Jfc — 1)!! 

2m- (2/ + 1) ~~ ~7r(2Z + 1)!! ^ 22 - 2fc + 1 (2fc)!! ' 

1 = 0,1,2,..., (4.14) 

which is proved in Appendix O 

Therefore, by the regularization of the divergent sum (j4.10|l . we exactly reproduce Hutson's 
solution (|4.4|) for x ~ for o(l). This coincidence is reasonable since we have started from the 
expansion of g(x) around x = 0. 

We may also calculate the physical quantities from (|4.11() . At our consideration, to the A 
order, we can integrate (|4.11(l for the whole interval [—1,1]. Indeed, the corrections come from 
near the edges, x <G [—1, — 1 + A] U [1 — A, 1], and g(x) is estimated to be of the order 1/VA there. 
Then, the corrections are of the order y/X, and thus can be neglected. 

The case of equally charged disks (upper signs in l|4.1|) and l|4.2fl ') has not attracted much 
interests. So far, no prescription has been proposed to solve 1)4. 2|l asymptotically. Our solu- 
tion (j2~5T|) gives 

C=±/VK.= ^ + -^log^ + 0«. (4.15) 

To obtain the last expression, the regularization of the divergent sum b in the same manner 
as <|4.10ll is required: 

OO OO _j 1/1 

b = Y —— = lim y ——k 2m = - log — , (4.16) 
^2™+! feVi^2m+l 2 A y J 

n—0 m—0 

where (|4.13f) has been used to obtain the last expression. The first term in (|4.15|l was also 
referred to in the Kirchhoff's paper. 25 The second term agrees with the analysis by Gaudin 7 
except 4 in the logarithmic function. 




(2m- 1)!! 
(2m)!! 



5 Concluding Remarks 

We have presented solutions of the Gaudin integral equation (|1.5|) . We have solved the integral 
equation in the form of power series. Based on the solution, the strong coupling case is analyzed, 
in relation to the Licb-Linigcr integral equation for the bosons with repulsive interaction. Then, 
we have a unified view on (^-function attractive fermions and repulsive bosons in one-dimension. 
Main results are summarized as follows. 

i) As to the thermodynamic quantities, the expressions in the strong coupling cases \"f\ 1 
are identical through the relations 1)3. 7|l . 

ii) The limit 1/7 = yields the Girardeau gas, 26 in the sense that the energy spectra coincide 
with that of free fermions. 

iii) Around the point 1/7 = 0, the system behaves like a hard-core bose gas with a diameter a 
given by (|3.11[) . This indicates the existence of dimers of tightly bound fermions ( 1 /-f = —0, 
a > 0), which can be regarded as bosons with strong repulsion (1/7 = +0, a < 0). 



15 



iv) The above regime of 7 can be smoothly connected to the weakly attractive fermi gas and 
the weakly repulsive bose gas at both ends with sufficiently small value of |^/| . This reflects 
the analyticity of the Love equation (|4.2fl in A on the whole real axis excluding the point 
A = 0. 

In the weak coupling case, we have again met the difficulty of divergences. Obviously, the 
solution 1|2.31[) holds only for x G [— 1 + A, 1 — A] (see Appendix [XJ. Near the edges \x\ w 1, a 
different method might be necessary to determine the precise behavior of the solution. Further, 
the relation between A and 7 again is important, since the regularization of A by 7 yields the 
finite expressions of the physical quantities. Up to now, our method is powerful enough to obtain 
the first two terms of the asymptotic expressions. We believe that other physical quantities, 
such as the dressed charge and the dressed energy, and other quantum integrable systems can 
be studied by use of this method. Those subjects remain as future problems. 



Appendices 



A Another Derivation of J22D and ffZTM. 



In this appendix we present the expansion of the integral, 

A f 1 . F(y) A f 1 



A f 

J(x, A) = — / dy ■ 

7T J-l 



dy 



m=0 



A 2 + (x - y) 2 7r y_ a A 2 + (x - y) 2 ' 
in power series of A for the case A <C 1. Since 

1 1 ^ (_i)«A 2 ™ 



1 



(x - yf 



A 2 + (x ~ yf 



= < 



A 



E 

n=0 



(x - y) 



2n+2 



1 

^ (-l)"(x - y) 2n 



E 



, 2 \2n+2 
X ~ V \ n=0 



, if \x — y\>\, 
if \x — y\<\, 



for a fixed x € [—1 + A, 1 — A], (|A.1J) becomes 



J(x,\) 



A 



■l,x-A)U(a;+A,l] J (x-\,x+X) 



dy 



A 2 + (x - y) 2 ^- 

v y ' m=0 



E 



a m y 



2 EEE a2 ' {2m ^ { ~ 1)namx2{m ~ l) 

m=0 n=0 1=0 
00 oc 2m 



21 J 21 - 2n - 1 

n „ (T ,2m-fc 



(A.l) 



(A.2) 



1 /9m\ f_1^n T 2™-1: 

^ E E E A2n+1 fc fc - 2 ; - 1 [ ( - 1)fe(1 + x)k ~ 2n " + (1 - s) " 2 "" 1 ] 



m=0 n=0 k=0 
oc m — 1 



rn—l n— 

„ oo oo m 

m=0 n=0 1=0 
1m 



— i a m a; 



2(m-n)-l 



log 



1 + X 

1-x 



l(m _ f 2m\ (-l)"a m a: 2 ' 

2Z / 2n + 2m - 21 + 1 ' 



(A.3) 



In the above, ^ means that we exclude the k = 2n + 1 terms in the sum with respect to k. 

k=0 



16 



Let us rewrite (|A.3|) so that we can easily see the dependences of A and x. For the purpose, 
we change the variables and the orders of the multiple sums. To begin with, we consider the 
first term of the r.h.s. of 1|A.3|) . The order of the triple series can be changed as 



oo oo m oo oo oo 



EEE = EEE- (a- 4 ) 

m=0 n=0 1=0 i=0 m=l n=0 

With a new variable p = m — I instead of m, the first term of the r.h.s. of (|A.3|) can be written 



as 



2 E E E x21 (2m ^ (- 1)namxKm ~ l) 



m=0 n=0 1=0 



tE^ e 

Z=0 lp=0 



2W 21- 2n - 1 



2Z + 2p\ 
21 



jai+p 



E 



(-i) n 



21 - 2n - 1 



z 2 p 



In exactly the same way, the fourth term becomes 



~ oo oc m 

!eee^ 2( 

m=0n=0 i=0 



m _ f2m\ (-l)"a m x 2 ' 

2Z / 2n + 2m - 2/ + 1 



p=0 l> 1=0 

The third term is rewritten as follows 

oc m — 1 



2p + 21 
2p 



a p+ i 



E 



(-iy 



2n + 2p+l 



m=l n=0 v 7 



= -_Va 2 " +1 (-1)' 

7T ^- — > 



n=0 



= -^E(- 1 )" A2n+1 iE 



oo oo 

E E 

_m— n+1 p— 

E 



2m 



1 



2(m-n+p) 



n=0 



.i=l 



2n + 1/ 2p + 1 

21 + 2n\ a,n+i 



2n + 1 J 2q - 21 + 1 



The second term is rather complicated. Since 



(-l) fc (l + x) k - 2n - 1 + (1 - xf- 271 - 1 = (-l) fc E t 1 + 



1=0 



r 2 9 



k - 2n - 1 



^E 



p=0 
fc/2-n-l 

2 E 

p=0 



2n - k + 2p 
2p 



r 2p 



jfe - 2n - 1 
2p 



-2/' 



^E 



2n-k + 2p+l\ 2p+1 
2p+l 

2 E 



p=0 



2p+ 1 



if < k < 2n and k is even, 



if 2n + 2 < fc < 2m and k is even, 



if 1 < k < 2n — 1 and k is odd, 



2p+1 , if 2n + 3 < & < 2m - 1 and k is odd, 



we decompose the sum with respect to k into four parts: 

2m 2n 2m n n— 1 



E' = E+ E=E+E+E + E 

k=0 k=0 k=2n+2 j=0 j=0 j=n+l j=n+l 

k=2j k=2j+l k=2j fc=2j+l 



(A.5) 



(A.6) 



(A.7) 



(A.8) 



(A.9) 
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For a fixed n > 0, the sum with respect to m can be decomposed into two parts: 

oo n oo 

E = E + E • 



(A.10) 



m— 771—0 m— n+1 

With and (|XlO)l . the double series with respect to m and k become 



oo 2m 



m m- 



-X \ 



m=0 fc=0 m=0 



-X \ 



EE' = E E + E + E E + E + E + E 



3=0 3=0 . 

\k=2j k=2j+l/ 



(A.ll) 



m— n+1 j— J— j— n+1 j— n+1 . 

\fc=2j fe=2j + l fc=2j k=2j + l J 



Now we proceed to evaluate the second term of the r.h.s. of (|A.3(I using (fO)l and l|A~TTJl : 



oo oo 2m / 9m \ / nn. ™2m-fc 

£ E E E A2n+1 ( T J fc _ 2 ; _ i K -1 )^ 1 + *) fc ~ 2 "" 1 + - ^~ 2ri ~ 1 ] 



m=0 n=0 fc=0 



k m— p=0 



/ 2m\ /2n - 2j + 2p 



V2j 



.3=0 

771 — 1 

+ E 

J=0 



2p 



-.r 



2m-2j+2p 



2j - 2n - 1 



2m \ / 2n - 2j + 2p\ 1 2 m-2j+2 P 
2j + l)\ 2p+l J2j-2n 



V 



E a ™E E 

m— n+1 p=0 Li— 

n-1 

+ E 

m n— 1 

E E 

77i—7i+l Lj— n+1 p— 

m— 1 j— n— 1 



2m\ (2n - 2j + 2p 



1 



2j -2n-\ 



2m-2j+2p 



2m ^(2n-2j + 2p^ 1_ 



2p + 1 ) 2j - 2n 



E 



2m\ (2j-2n-\ 



2j 



2p 



1 



2m-2j+2p 



2j - 2n - 1 



/ 277i \/2j-2rc\ 1 2m - 27 - +2p 



-y(-i) n A 2n+i 



n=0 



J =71+1 p— 

oo n 2g 

EE E 

Lr=0 9=0 2=0 
n — 1 n 2q 

E E E «*(-d 

r=0 q=r+l l=2r+l 



2q\ (2q - 2n - 1 - I 
2r-l 

2q\ (2q-2n-\-\ 

2r-l 



1 



2q - 2n -I - I 
1 

2g-2n- 1 -V 



OO 



2r + 2n\ 1 



'-\-n 

7-1 

oo r — 2 2 r 



2r y -1 



..2r 



+EE E a 9+»+i(- 1 ) 

r=2 q=0 Z=2g+2 

oo r — n—2 2r 

- E E E a ?+n+i(- 1 ) 

r=n+2 q=0 /=2<?+2n+3 
oo oo 2r 



2g + 277 + 2\ /2g + 1 - / 



1 



2r - I J 2q + 1 - I 



2q + 2n + 2\ (2q + 1 - I 



r=0 q=r Z=0 



2<? + 2n + 2\ (2q + 1 - I 



I 



2r 
1 



,2r 



1 



2q+l-l 



2r - I J 2q+l-l 



(A.12) 



To obtain the last expression from the middle, we have again changed the variables and the 
orders of the multiple sums. For instance, in the hrst square bracket [ ], we have replaced m, j 
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and p with q = m, I = m — j, r = m — j + p, and used the relation 
n + k — 1 \ , _ Nfc ( — n 



(-1) "I k j for n>l, fc>0, (A.13) 

to collect the terms of 21 and 2/ — 1. 

The combination of (|A.5jl . (|A.12jl . (|A.7jl and ljA.6|) yields the following form of (|A.3jl : 

r) °o r> oo 



J(x, A) =- V X 2n [A n (x) + B n {x)} + - V(-l) n A 2n+1 [C n (x) + D n (x) + E n (x) + F n (x)] , 

7T ' — ' 7T » — ' 

n=0 n=0 

(A.14) 

where A n (a;) , £?„ (x) , C„ (ar) , D n (x) , (x) and F„ (x) are defined as 

M*>--±±(\ + *)£^-^, (A,5, 



p__0 g=0 

oo oo 



p=0 9=0 

n / p 2q n 2p \ oo n 2q 



2p 



npqr^ j 



p__0 \q— r— g— p+1 r— / p— n+1 g— r— 

(A.17) 

oo / p— 1 2g n+p— 1 2p \ oo 

^ ^ I ^ ^ ^ ^ ^q-^npgr*^ ^ ~l~ ^ ^ ^ ^ G-g-^npgr^ ^ J ~r~ ^ ^ ^n-}-p^np.p-^n,2p^ ^; 
p— 2 \g— n+1 r—lq — ln, Q—P r— 2g — 2n / p— 1 

(A.18) 

oo oo 2p 

E n{x) = Y J £ J]MwX 2p , (A.19) 

p— g— n+p-|-l r— 

Here, the sequence 

/2g\ /2g - 2n - 1 - r 



A ^^(-lY\ r )[ 2p _ r ) 2q _ 2n _ l _ r > n, P , q ,r & Z> , (A.21) 

has been introduced to simplify the expressions. Note that sums in the above whose lower 

p 

bounds are greater than its upper bounds, such as ^ , are treated as zero. 
From (|A.14|I . we obtain the expansion of J(x, A) with respect to A as follows: 

oo 

J(x,\) = Y / * n Mx), (A.22) 

n=0 

where the first few J n are 

Jo(x) = -[Mx)+B (x)], (A.23) 

Ji(x) = - [C (x) + D (x) + E (x) + Fq(x)} , (A.24) 

J 2 (x) = -[A 1 (x)+B 1 (x)], (A.25) 

J 3 (x) = --[C 1 (x)+D 1 (x)+E l (x)+F l (x)]. (A.26) 
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From 1A.23I). llA~25l) and 



E 

n=0 



(-1)" = 7T 

2n+l 4 : 



we easily verify that 



J (x) =Y,a p x 2p = F(x), 

p=0 



(A.27) 

(A.28) 
(A.29) 



Mx) = ~ + l)(2p + l)a p+1 x 2 *\ 

p=0 

Use of l|A.28(l in the Gaudin integral equation 1)2.3(1 leads to ((2.24(1 . To simplify the expressions 
of J\(x) and 1/3 (x), we need a little calculation and some formulae. The results are, 



7T * 



E 



2p+ 1 



p=0 



2n - (2p + l) n 
/2p + 3 



p=0 Ln=0 

00 



^ 2n -(2p + 3) 



r-2P 



(A.30) 



(A.31) 



Below, we shall explain these shortly. We observe that equation ((A. 30(1 is in accordance 
with (|2~2S|) . 

Therefore, we conclude that our formal expansion ((2.7(1 is also valid for small A in the interval 
[-1 + A.1-A]. 



Calculation of J\ (x) 



From ((A. 24(1 with l(A.17() - l(A.21(l . we have 

00 00 p 00 00 2p 

-zJi (x) = ^2 a A 0p oox 2p + ^2 E a q-Ao pq ;2qX 2p + E E E a q-Ao pqr x 2p 



p=0 



p=l q=l 



p—0 q—p+1 r— 



00 p 



2q\ a 



q - x 2p 



yy , 

OO ^ OO p ^ 



,2p 



p=0 



00 00 



p=l 9= l 
2? 



E E a ? 

p=0 g=p+l 

00 p 

EE 

p=l g=l 

00 p 



E( ^ ( r )( 2p-r ) 2<? 1 



r 2 '' 



2q 



2q-2p-l 



,.2p 



EE— a ^ 2p +E E 2q -(2p+i 

p=0 9=0 p=0g=p+l v ' 



00 p 



n^ 2p +EE 



29 



p tl^2g-(2p+l) 



EE 



2p+ 1 



2« -(2p + !)"«- 



(A.32) 



p=0 g=0 

Here, to obtain the third expression from the second, the following formula has been used: 



EM) 



1=0 



m + 1\ (m — l\ 1 



/ J \p — I J m — I 



(-1) J 



p+l\ 1 



1 J m — p 



m > p + 1 . 



(A.33) 
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Proof of 
Set 



P P (m) 
TO — p 



i=0 



m + 1\ /to — Z\ 1 



Z J \p — I J m — I 
It can be seen that 

P p {m) is a polynomial of to whose degree is equal to or less than p. 
It will be shown that P p (m) is a constant, and actually, 

P p (m) = (-iy(p+l). 



(A.34) 



(A.35) 



(A.36) 



When p = or p = 1, we can verify it directly, so we assume that p > 2. 

p 

First, let us denote the Z-th term of the sum ^ in the r.h.s. of (IA.34|) as Qi(m)/(m — p), 

1=0 

namely, 

p 

P p (m) = ^Qi{m). (A.37) 

1=0 

Then, Qi(m) is explicitly written as 



Qi(m) = (-1) 



TO -(- l\ / TO — Z\ TO — p 



Z J \p — I J m — I 
' (to — 1)(to — 2) • ■ ■ (m — p) 

( 1 y ( m + !) m ■ ■ ■ (™ ~ Z + 1)(to - Z) ■ ■ ■ (to - p) ifl< ; <p _ 1 



if Z = 0, 



?!(p-0! 

( 1 y P ( m + 1 ) m '"( m ~-P + 2 ) 



(A.38) 



if I = p. 



Here, • • • ( • ) • • • denotes that such term is excluded from the sequence of the product. 
From (TOTl) and (lA~38l) . P p (0), P p (l), . . . P p (p) are calculated as follows: 



P p (fc) 



Q fc (fc)+Q fc+1 (fc), if0<fc<p-l, 



\Qp(p), 
(-i) p (p+i). 



if fc = p, 



(A.39) 



Therefore, (fQ9)l with IpQgjl assures that P p {m) = (— l) p (p + 1) for all to, which completes 
the proof. 



Calculation of Jz{x) 

This could be done in almost the same, but more complicated way than JTi(x). We sketch 



21 



the calculation. From ifp6)i with JS33-(S23J, 

wc have 



oo 2 



p=0 

oo p 2q 



p— 1 r— 



+ QqAlpqrX P 

p=2 q=2 r=2q-2 

OO P+l /r) s 



~ /2p + 2\ 1 



p=i 



p=lg=2 

oo oo 



tS. x 2p 



3 J 2q - (2p + 3) 



E E a </ 

p=0 g=p+2 



2p 



a oE 

p=0 
oo p+l 

+EE 
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In the above, the second expression is obtained from the first by applying the identity 



EM) 



(=0 



m + 3\ fm — l\ I 



I J \p — I J m — I 



-1) P ( P+ J)^-, m>p + l, 



3 / m — p 



(A.40) 



(A.4I) 



for the sum with respect to r in the last term. We note that (|A.41|I can be proved in a similar 
way as (|A.33I> . 



B Solutions of the Lieb-Liniger and the Gaudin integral 
equations around 1/A = 

In this appendix, we solve an integral equation, 

nl H(y,r])dy _ 1 



H{x,rj) — — f - 

7T J -I 1 



+ i1 2 (x~ y) 2 2ir' 



(B.I) 



around 7] = 0. We remark that, H{x,rj) gives the solutions F(x) of (|2.3(l and g(x) of (|3.1(l in 
the strong coupling case (i.e. |1/A| <C 1), 



F(x) =4wH(x i -l/X), 
g(x)=H(x,l/X). 



The physical quantities are expressed as follows, 
for 7] < 0, 



7f 



4?7 J H(x,r])dx, 



Q = D F 



H(x, rj) dx 



3 A 



1 ~ 3 

H(x,rj)dx 

-l 



(B.2) 
(B.3) 



(B.4a) 
(B.4b) 
(B.4c) 
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for r\ > 0, 



7b 



7] J H(x, rj) dx 



K = D B J H(x, i]) dx 
^ = D'i [ x 2 H(x, rj) dx 



1 I -3 

H(x, if) dx 



(B.5a) 
(B.5b) 
(B.5c) 



Hereafter, we denote H(x,rj) as H{x) for simplicity. With the expansions of H(x) and the 
kernel. 



H(x)=Y. b P x2P > 



p=0 



1 



equation l|B.l() becomes 



1 + if(x - yf 



oo oo oo 
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E^-±£££(-ir 



2n+l "P 
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r 2(n-0 _ J_ 
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which, as in §12 reduces to the following set of algebraic equations for {bp}: 
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From 1B.9|1 . we get 

bo = 

h = 



1 

2T 



2 4 2 2 / 4 1 



7T 7T^ 7T \ 7T 



2/20 1 16 



1 + -TJ-2 1 



7T V 371" 2 5 7T 4 



0(r? 6 ), 



& 2 = + 0(»? 6 ), 



& p = 0(^ p+1 ), p>3. 
By use of the above solution, 77 may be expressed as 



where 



In terms of £, we have 



1 



7b + 2 4 7F + 2 



1 



2tt(1 - 2£) 



i-2€ 

6 2 = 7r 3 e 5 + o(? 6 



1- +o(£ 6 ), 



o(£ 6 ), 



(B.6) 
(B.7) 

(B.8) 

(B.9) 



(B.10) 

(B.ll) 

(B.12) 
(B.13) 

(B.14) 
(B.15) 

(B.16a) 

(B.16b) 
(B.16c) 
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It follows from the above relations that for £ > 0, 



K = irD B (l-2$ ( 1 + -7r 3 e 



O(f). 



(B.17) 
(B.18) 



Other physical quantities are readily obtained, which reproduce . 5|) and (|H.6|I . 

C Proof of (BHD 

We begin with the following expression of the complete elliptic integral of the first kind, 



K{k) 



d0 



IE 



'o \J\ - k 2 sin 2 9 2 m=0 
If we consider the expansion such that 

2 r 1 dt n d0 
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(2m - 1) 



(2m)! 
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o yl - t 2 sir 



V] c n i", for x > 0, 



(C.l) 



(C.2) 



n— — oo 



we see from (|C.1|) that 



Co 



E 



(2m - 1) 
(2m)!! 



1 



2m - 21 - 1 



(C.3) 



In the l.h.s. of l|C.2jl . since the double definite integral is convergent for x > 0, we can invert 
the order of the integrations. To integrate with respect to t, we change the variable from t to 
ip, where tp is defined as <p = arcsin(t sin (9) , < <p < tt/2. Then the integration in the l.h.s. 
of (|C.2|I becomes 



2 f 1 dt 

IT 



t 2l+2 



dO 



o s/l — t 2 sir 



2 n 

7T 
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^ arcsm 
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By use of the formula, 



(C.4) 
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{2l)\\ cos<p y> (21 -2k- 1)!! . 
~(2Z + l)!!sin 2/ +V^ (2Z-2A:) 



II 



sin 



(C.5) 



the r.h.s. of 1C.4I) is rewritten as 



2 PrfJ » 



cos 



E 

fc=0 



(2Z - 2fc- 1)!! 2fe 
(2I-2A)!! 



x 2 sin 2 



E 



fc=0 



,(«-2*-l)H^ M( 



(2/ - 2jfe)!! 



(C.6) 



Since the second term in the square bracket [ ] consists of odd powers of x, only the first term 
contributes to cq, which gives 



Co 



2 (20!! (2i - 2k - 1)!! 1 
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(C.7) 
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Therefore, we obtain 



OO r 



(2m- 1)!! 
(2m)!! 



1 



2 (2Q!! 
tt(2Z + 1)!! 



(21 -2k- 1)!! 1 
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2m - (2/ + 1) 



E 



2 (21)!! 
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l 



1 (2fc-l)!! 
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E 



Z = 0,1,2, .. 



(C.8) 



which is nothing but Q4.14p . 
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